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Abstract. 

The classical problem of irrotational long waves on the surface of a shallow layer of an ideal fluid 
moving under the influence of gravity as well as surface tension is considered. A systematic proce- 
dure for deriving an equation for surface elevation for a prescribed relation between the orders of 
the two expansion parameters, the amplitude parameter a and the long wavelength (or shallowness) 
parameter /?, is developed. Unlike the heuristic approaches found in the literature, when modifica- 
tions are made in the equation for surface elevation itself, the procedure starts from the consistently 
truncated asymptotic expansions for unidirectional waves, a counterpart of the Boussinesq system of 
equations for the surface elevation and the bottom velocity, from which the leading order and higher 
order equations for the surface elevation can be obtained by iterations. The relations between the 
orders of the two small parameters are taken in the form fj = 0(a n ) and a = 0((3 m ) with n and m 
specified to some important particular cases. The analysis shows, in particular, that some evolution 
equations, proposed before as model equations in other physical contexts (like the Gardner equation, 
the modified KdV equation, and the so-called 5th-order KdV equation), can emerge as the leading 
order equations in the asymptotic expansion for the unidirectional water waves, on equal footing 
with the KdV equation. The results related to the higher orders of approximation provide a set 
of consistent higher order model equations for unidirectional water waves which replace the KdV 
equation with higher-order corrections in the case of non-standard ordering when the parameters a 
and /? are not of the same order of magnitude. The shortcomings of certain models used in the lit- 
erature become apparent as a result of the subsequent analysis. It is also shown that various model 
equations obtained by assuming a prescribed relation (3 — 0(a n ) between the orders of the two 
small parameters can be equivalently treated as obtained by applying transformations of variables 
which scale out the parameter /? in favor of a. It allows us to consider the nonlinearity-dispersion 
balance, epitomized by the soliton equations, as existing for any /3, provided that a — > 0, but leads 
to a prescription, in asymptotic terms, of the region of time and space where the equations are valid 
and so the corresponding dynamics is expected to occur. 
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1. Introduction 

The behavior of surface gravity waves on shallow water has been a subject of intense 
research. In particular, the famous Korteweg-de Vries (KdV) equation, which is the 
prototypical example of an exactly solvable soliton equation, was first introduced as 
a unidirectional nonlinear wave equation obtained via asymptotic expansion around 
simple wave motion of the Euler equations for shallow water. 

The system of equations describing the long, small-amplitude wave motion in 
shallow water with a free surface PQ [I] involves two independent small parameters: 
a, which measures the ratio of wave amplitude to undisturbed fluid depth, and (3, which 
measures the square of the ratio of fluid depth to wave length, and no relationship 
between orders of magnitude of a and (3 follows from the statement of the problem. 
The KdV equation 

3 1 

Vt + Vx + 2 ai lVx + qPV3x = (1) 

emerges at first order (in both parameters a and 0) in the asymptotic expansion as 
an equation for the surface elevation 77 associated with the right-moving wave. The 
derivation assumes (sometimes tacitly) that /3 = 0(a). It is evident that in the case, 
when a and /3 differ in their orders of magnitude, the leading order equation maintaining 
the balance between linear dispersion and nonlinear steepening, which is the primary 
physical mechanism for the propagation of solitary shallow water waves, should change 
its form. The same holds true for the equations which (like the higher-order KdV 
equations) address higher order effects. 

A heuristic approach to deriving model equations for unidirectional water waves is 
frequently used when some additional terms are included into the equation for the surface 
elevation based on relations between the orders of parameters. However, this may lead to 
inconsistencies. For example, the assumption a > (3 > a 2 is made in [5] while the terms 
involving a, (3, a 2 , a(3 and (3 2 are kept and, accordingly, the terms involving a 3 , a 2 (3, af3 2 
and (3 3 are neglected. It is readily seen that the relation (3 = a r with 3/2 < r < 2 which 
satisfies the above inequality is in conflict with the truncation made: the neglected term 
~ a 3 is as important as the retained term ~ f3 2 . In [6], based on an inequality 0(f3) < 
0(a), the truncation is made such that the terms involving a, (3, a 2 , a(3, a 3 and a 2 j3 are 
kept and the terms involving a A , a 3 (3 and f3 2 are neglected. However, this choice of trun- 
cation is questionable since there exists no relationship between the orders of (3 and a 
of the form (3 = 0(a r ) (or a = 0(f3 r )) for which such a truncation is consistent. Indeed, 
assuming (3 = 0(a r ) with r > 1 (which is compatible with 0(f3) < 0(a)), one can see 
that the two requirements a 4 < a 2 (3 and (3 2 < a 2 (3 lead to conflicting results: r < 2 and 
r > 2 respectively. Thus, such a heuristic approach does not provide a reliable way to 
determine even a form of the equation for surface elevation and, what is more, it does not 
allow determining coefficients of the equation. It is well known that the solution prop- 
erties may strongly depend on the relations between the coefficients - the higher order 
KdV equations can be mentioned in this respect (see e.g. [I] [9] and references therein). 
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In general, to arrive at a consistent model equation for water waves, the ordering 
of terms should be made in the original asymptotic expansion for unidirectional water 
waves based on a prescribed relationship between orders of magnitude of a and (3. 
Then a consistent truncation of the expansion can be made and the related leading 
order and higher order evolution equations can be defined. In the present paper such a 
procedure for deriving an equation for surface elevation for a prescribed relation between 
the orders of the two expansion parameters a and (3 is developed. It makes possible a 
systematic study of different particular cases and corresponding leading order and higher 
order equations. The following special cases are considered: (3 = 0(a 2 ), a = 0((3 2 ), 
(3 = 0(a 3 ) and a = 0((3 3 ). The analysis is aimed at deriving an equation for the 
surface elevation having a form of an evolution equation; therefore, equations which, 
like the Benjamin-Bona-Mahoney (BBM) equation [10J , contain the time derivatives 
in the higher order terms are excluded from consideration. The results of the analysis 
show, in particular, that some evolution equations proposed before as model equations 
in other physical contexts can play the role of a model equation at the leading order of 
the asymptotic expansion for the unidirectional water waves on equal footing with the 
KdV equation. Some of these equations, both integrable and non-integrable, are known 
to have a rich structure of solitary wave solutions which differ in their properties from 
the KdV solitons. Thus, the leading order soliton dynamics in the unidirectional water 
wave problem can differ from the one described by the KdV equation. The equations 
obtained in the higher orders of approximation, in general, also differ from the KdV 
equation with higher order corrections. It is worth noticing that the above differences 
from the standard model are not due to taking the surface tension into account. New 
equations and dynamics arise even in the classical formulation, when capillary effects 
are neglected, if the ordering is non-standard ((3 and a are not of the same order of 
magnitude). Including surface tension in general does not alter the structure of the 
leading order and higher order equations, only some specific cases, like the case r = 1/3 
of the standard analysis, should be considered separately. 

The paper is organized as follows. In Section 2 following the Introduction, we 
present the statement of the problem, the basic equations and the outline of the 
procedure. The main ideas of the analysis are described in more detail in Section 3 
where the procedure is presented for the best studied case of (3 = 0(a). The cases when 
the relation (3 = 0(a) does not hold are studied in the subsequent Section 4. In the said 
section the analysis is restricted to the pure gravity waves in order to better explain the 
main points and also to demonstrate that the differences from the standard model are 
not due to taking the surface tension into account. The results for the gravity-capillary 
waves are listed in the Appendix. In Section 5, the concluding remarks are given and an 
alternative interpretation of the results, based on a transformation of variables which 
scales out the parameter (3, in favor of a, is discussed. 
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2. Outline of the procedure 

Consider the standard system of equations describing the two-dimensional irrotational 
wave motion of an inviscid incompressible fluid in a channel with the flat horizontal rigid 
bottom and the free surface under the influence of gravity as well as surface tension. 
After an appropriate choice of non-dimensional variables, the equations of motion and 
boundary conditions can be reduced to the system written in terms of the velocity 
potential <j)(x,y,t) and the surface elevation rj(x,t), see e.g. pQ: 

P<p xx + K = 0, < y < 1 + a V (2) 
<P y = 0, y = (3) 
T) t + a(f) x r) x - —<j> y = 0, y = 1 + ar] (4) 

1 1 OL Tj 

& + 2 a<Pl + 2~/y + v ~ T \i + Jp V 2f* = °» y = l + ^ ( 5 ) 

where t is time and x, y are respectively horizontal and vertical coordinates, with y = 
being the bottom. The non-dimensional variables are defined as follows (after non- 
dimensionalizing, the tildes have been omitted): 

x= v y= w ^=a' t= i7^' 0= wmvm ( ) 

where g is the acceleration due to gravity, H is the upstream mean depth and a and 
L are typical values of the amplitude and of the wavelength of the waves. Equations 
(EJ) — (EI) contain three non-dimensional parameters: the amplitude parameter a — j|, the 
wavelength parameter (3 = and the Bond number r = -^p, where T is the surface 
tension coefficient and p is the density of water. 

Equations and are satisfied by making a standard substitution 

^ = S (2m)! <9x 2m ^ ' 

where f(x,t) = 4>\ y =o- Substituting (J^J) into the surface conditions (BJ and (EJ) and 
differentiating ([5]) with respect to x yields a system of two equations for the surface 
elevation T)(x,t) and the horizontal velocity at the bottom w(x,t) = f x in the form of 
infinite series with respect to /3. We are interested in considering weakly nonlinear small 
amplitude waves in a shallow water, so we will treat a and /3 as small parameters. 

In the zero order in both a and /3, the system of equations for w and r\ reads 
Vt + w x — 0, w t + r] x = so both w and rj satisfy the linear wave equation 
Ctt — Cxx = which describes waves traveling in two directions. A wave moving to 
the right corresponds in this order of approximation to w = n and rj t + r] x = 0. To 
derive the equations describing right-moving waves in higher orders in a and 0, we can, 
along the lines of PQ, reduce the system of equations for w and r] to an asymptotically 
equivalent set of equations consisting of a relationship between the horizontal velocity 
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w and the surface elevation 77 and an evolution equation for the elevation. To do this, 
we set 

oo 

w=J2 Rij^ft, (8) 

i,j=0 

where Rij depend on r\ and its ^-derivatives, and possibly some nonlocal variables, with 
-Roo = Vi an d require that rj satisfy an evolution equation of the form 

oo 
i,j=0 

where Sqq = —r] x and in general depend on 77 and its x-derivatives. The functions 
Rij and Sij are determined from the requirement of consistency of (jHJ) and ([9]) with the 
above system of PDEs for w and rj. To implement this, an iterative procedure starting 
from the zero order of approximation and continuing to the higher orders is applied, see 
the subsequent sections for details. In each order, the t-derivatives of 77 are replaced by 
their expressions from the lower order equations. 

However, it is obvious that truncating the asymptotic expansions and keeping only 
the terms up to certain order requires the knowledge of relationship between the orders 
of magnitude of the two small parameters a and ft, because otherwise it is impossible 
to determine which terms should be retained and which can be neglected. A commonly 
used assumption is that a and ft have the same order of magnitude (ft = 0(a)). Then, 
choosing, for example, a to be a primary parameter and retaining the terms up to 
0(a n ), we arrive at the so-called n-th order Boussinesq system p], If the first order 
Boussinesq system is considered and the corresponding order expansions are taken for 
dHJ) and (jHJ), then equation (jUJ) for the elevation 77 takes the form of the KdV equation 
dlj, cf. pQ, see also the next section. The same procedure continued to next orders 
results in the KdV equation with higher order corrections. If the relationship ft = 0(a) 
does not hold, then one needs an alternative assumption relating the orders of a and ft 
to make a truncation of the expansions consistent. 



3. Procedure for the case of ft — 0(a) 

In order to explain how the forms of the expansions (jHJ) and are determined up to a 
certain order through an iterative procedure, we will first present the procedure for the 
best studied case of ft = 0(a) (we can set ft = a without loss of generality). We will 
also consider the problem without surface tension to make the analysis as transparent 
as possible. Then equation (j5J) is replaced by the following 
1 la 

& + -a4>l + --p4>l + n = Q ) y = l + a V . (10) 

If in the system of equations for w and 77 the terms in the second power of a 
are retained and the higher order terms are dropped, we arrive at the second order 



Ordering of parameters in the water wave problem 



6 



Boussinesq system 

r] t + w x + a ((iwi) x - ^V3x^j + « 2 y~ (v™2x) x + ^) W ^ = °> ( n ) 



Vx + a (ww x - ^2xt \ + a 2 (~ (v w xt) x + ^w x w 2x - ^ww 3x + i^Wzxt) = 0. (12) 

In the lowest (zero) order, the system (fTTjl . (1T21 reads r\ t + w x = 0, w t + rj x = 0, 
and the equivalent system (jSJ), (j5J) describing a right-moving wave is reduced to 
w = f], rjt + 7/2 = 0. In the next order iteration, we look for a solution for w corrected 
to first order as 

w = r] + aQ {1) , (13) 

where is a function of r] and its x-derivatives, and substitute (1T3]) into equations 
(TTT|) and (IT21) with the terms of order higher than 0(a) dropped. Upon the substitution, 
the equations in question become 

Vt + Vx + a(2rjr] x - -T] 3x + Q«) = 0, (14) 

Vt + Vx + a(2r]T] x - ~r) 2x t + Qt 1J ) = 0. (15) 

The function is sought such that the two equations CHI) and (IT51) agree (up to the 
first order in a) upon expressing all the t-derivatives of rj in terms of its x-derivatives 
using the zero order equation r\ t + r\ x = 0. This yields = —\rf + ^r] 2x . Then 
equations (TT3"|) and (ITB1) become 

iu = rj + «(-^ 2 + 3*723) > r)t + Vx + ®[-r)r) x + -r) 3x ^ = 0. (16) 
The equation for 77 is reduced to the KdV equation in a standard form 

Vi + SWx + V3x = (17) 
by the change of variables (x,t) — > (x,t), where 

£= \ZI (x_t) ' £= ivf ai (i8) 

At the next step, the above expression for w is corrected to second order in the 

form 

w = rj + a(--ri 2 + -r^ + a 2 Q {2) . (19) 

It is substituted into (TTTT) and (fT2l . and then all the t-derivatives of 77 are replaced by 

their expressions through the x derivatives using the lower order equation, namely, the 

second equation of (fl6j) . The condition of consistency of the two equations obtained in 

such a way leads to an equation for Q^ 2 \ a solution of which is expressed in terms of r\ 

and its x derivatives. As a result, we have 

/ 1 2 1 \ 2/I s 3 2 1 1 \ 

w = 77 + a (--77 + -r/ 2a; j+a ^ + — ^ + ^PTac + Jg^J, ( 20 ) 

/3 1 \ 2 ^ 3 2 23 5 19 \ 

Vt + Vx + ay-rjrj x + -rj 3x J + a (--77 77 x + — + —nn 3x + ^V5x J = 0. (21) 
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The equation for 7] can be reduced to the KdV equation with the first order correction 
in standard form. The procedure for determining the expansion (Q of the evolution 
equation for the right-moving wave can be continued to any order and yields the KdV 
equation with higher order corrections. 

The next (third) order corrections with the terms up to the seventh-order spatial 
derivatives included are given in [12], [X3J . The second and third order corrections (up 
to seventh-order derivatives) for the case of nonzero surface tension (r ^ 0) can be 
found in [H] (at the end of Appendix A) as a particular case of more general equations 
for bi-directional waves. Note that in sections 4 (for r = 0) and 5.3 (for r ^ 0) of the 
present paper corrections including the terms up to the ninth-order spatial derivatives 
are calculated for the case of a = 0(f3 2 ). 

It should be emphasized once again that our procedure is aimed at deriving 
equations for rj which have the form of an evolution equation, and, accordingly, all the 
i-derivatives in the terms of the order higher than zero are replaced by their expressions 
through the x- derivatives. Therefore, applying this procedure cannot yield equations 
which, like the BBM equation [10], contain the time derivatives in the higher order terms. 

It is also worth noticing that there exists a possibility to introduce certain freedom 
into the Boussinesq system. For example, a class of Boussinesq systems which are 
formally equivalent to the system displayed in (TTTj) (TT2|) can be derived using other 
variables instead of the horizontal velocity at the bottom w and employing the lower 
order equations in higher order terms [I], [15], [llj . It might seem that this freedom, 
revealing itself as free parameters present in the Boussinesq equations, should result in a 
freedom in the equation for the surface elevation derived from the Boussinesq equations 
under the assumption of unidirectionality. However, this is not the case: it can be readily 
checked that all those different but asymptotically equivalent systems are reduced to 
the same high order KdV equation for the surface elevation (1211) if the wave moving to 
the right is specialized. 



4. Examples of ordering 

In this section, examples of a non-standard ordering when the relation = 0(a) does 
not hold are studied. Two important special cases, (3 = 0(a 2 ) and a = 0((3 2 ), are 
considered in more detail. The analysis is restricted to the pure gravity waves in order 
to better explain the main points and demonstrate that the differences from the standard 
model are not due to taking the surface tension into account. The results for the problem 
including surface tension are presented in the Appendix. 

Starting from the case of (3 = 0(a 2 ), we first write down the system of equations for 
w and rj obtained by keeping all the terms of the order not higher than f3 2 , /3a 2 and a 4 : 

1 1 

r) t + w x + a {r)W) x - ~(3w 3x - -a(3 (r]w 2x ) x 

- l -a 2 p {v 2 w 2x ) x + = 0, (22) 
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w t + r]x + aww x - -(3w 2x t + a(3 ( - {nw xt ) x + -w x w 2x - -ww 3 



+ a 2 (3 ( w x (rjw x ) x - - {v 2 w xt ) x - w {j]w 2x ), 



o. 



(23) 



Next, we apply the iterative procedure described in the previous section to determine 
the form of the unidirectional wave equations (JS} and flU]) for the case of (3 = 0(a 2 ). 
The resulting equations, with the terms up to 0(a 4 ) retained, read (we have used the 
square brackets to gather the terms having the same order of magnitude) 



w - 

+ 

Vt - 
+ 



n 



a— + 
4 



+ 



,5rj 



a 6 — '- + a(3 — + 



V 2 V2x 3wil 3z 
8 32 16 

8 6 



64 



+ 



+ 



%nl , VV2* 



16 



2 Vix 

10 



a 3 V 3 Vx + a(3 ( ^VxV2x + t^Wsx 



rf x dx. 
/23 



(24) 



16 



15 



128 



4 4 i 2 a 

a T) r) x + a p 



23 



19 



Y^V V3x + —VVxVxx + 3 ,- } 



V24 

19 o2 
360 



12 



0. 



To make the things even more clear, rewrite the last equation taking 
(B = 0(1)) and ordering the terms according to powers of a. We obtain 



(25) 
Ba 2 



Vt + ilx + 7; ar 1Vx + a [ --T) n x + -Brj 3x 



+ a q 



2 

15 



+ a 6 



16 



128 



23 



19 



V 3 Vx + B 



19 



23 5 

^VxV2x + J^VV3o 



V Vx + B ( — fj rj 3x + —ririxrixx + t^Vx ) + ^ B rj 5x 



0. 



j^'' 'I'' 1 j J-j 'I'll'll } •'1(':r\^ ' / ' >r ' (^^) 

It is immediate that an equation involving both nonlinearity and dispersion is 
obtained at the leading order, which is now second in a and first in /3. Therefore 
this leading order equation contains an extra term — \a 2 rfrj x and reads as follows: 



Vt + Vx + -jOO]r\ x 



^a 2 r] 2 r] x + -/3r] 3x 



0. 



(27) 



Thus, if — 0(a 2 ), then the leading order equation is not the KdV equation but 
the Gardner equation which is a linear combination of the KdV and of the modified 
KdV equation. The Gardner equation has appeared in the literature in other physical 
contexts; in particular, it was derived in an asymptotic theory for internal waves in a 
two-layer liquid with a density jump at the interface [IE] , [T7] . Our derivation shows that 
the Gardner equation emerges in the classical water wave problem as the leading order 
equation in the case of (3 = 0(a 2 ). The Gardner equation is integrable and possesses 
solitary wave solutions but the Gardner solitons may differ in their properties from their 
KdV counterparts, see e.g. [T8] . 

The Gardner equation ([27)1 can be transformed into the modified KdV equation 

fk = fkz + Gfffk, (28) 

where fj is a shifted variable and (x, t) are the rescaled variables in a moving frame. 
Equation (J28"j) is well known to be integrable, see e.g. [21 Ej. In addition to standard 
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soliton solutions, it has solutions in the form of 'breather solitons' and also solutions 
describing breather-soliton interactions. In view of the fact that the transformation from 
( 1271) to ( 128]) includes a shift of the dependent variable, soliton solutions of equation (128]) 
for fj can be relevant for the original problem in terms of 77 if the flows with hydraulic 
jumps are considered. 

The form of the higher order corrections to the leading order Gardner equation is 
also evident from equation ( 125|) (or (126]) ). Note that equation (125]) has the differential 
structure of a combination of the Gardner equation ( |27j) and its first commuting flow; 
this feature is similar to what is observed for the KdV with a higher order correction in 
the case of (3 = 0(a). 

In P|, the so-called second and third order approximations of water wave equations 
are studied for the case 0(f3) < 0(a) of [19] specified to (3 ~ a 2 . The comparison of these 
equations with ( 125]) and ( 126]) shows that the ordering (and hence the truncation) used 
in [6] are invalid. In particular, in the second-order approximation equation the terms 
involving a/3 are present but the same order term involving a 3 is missing. Likewise, in 
the third-order approximation equation the terms involving a 2 (3 are retained but the 
same order terms involving a 4 and j3 2 are omitted. 

Consider now the case a = 0((3 2 ). Then the basic system of equations for w and 77 
obtained by keeping the terms up to 0((3 A ) (or 0(a 2 )) has the form 

\- ^a(3 2 (nw ix ) x + ^1^0 V» = (29) 



1 11 1 

Vt+w x - -(3w 3x + a (i]w) x + y^ 2 ^ - 2 (v™2x) x - 



24 

w t + Vx 
1 



1 1 

-f3w 2xt + aww x + — f3 2 w ixt + a(3 



+ 



720 
1 



P w&xt + ol^ 



7 (VW 3 xt) 

6 



E + —w 2x w 3x 



- (vw x t) x + t^w x w 2x 

1 1 

-w x w Ax + —ww 5x 
o Z4 



-ww 3x 



40320 



(3 4 w 8 . 



xt 







(30) 



An equivalent system of the unidirectional wave equations (|8]) and (|9]) truncated to keep 
terms up to 0(/3 4 ) is (the meaning of the square brackets is the same as in (125]) ): 



w 



+ 



1 1 



+ 



1 n / 9 x 61 

16 ' ' ; 1890 



P 3 V6x 



a 2 r/ 3 + a(3 2 



8 



1 



163 
360 



1091 7 



+ 



1261 
113400 



P A V8x 



Vt + Vx + ^PV3x + 

6 



3 19 02 

2 aVV * + 360^ ^ 



+ 



.23 5 

a P \ ^VxV2x + y^vvsx 



+ 



55 
3024 



+ 



^a 2 r] 2 r] x + a(3' 2 I ^-r] 2x r] 3x + 



317 



1079 



19 



288 



1440" ^ + 80™* 



11813 
1814400 



/3 4 V9x 



P%x 



0. 



It is immediate that for a 
nonlinearity and dispersion is 

1 3 19 

Vt + Vx + qPVsx + 2 ar lVx + 



0((3 2 ) the equation including at leading order both 



(3\x 



(31) 
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By the change of variables 





equation (!3Tj) can be reduced to the following 



19 

rji + Qririx + n 3x + Mn 5x = 0, M — —a. (33) 

This equation, which is frequently referred to as the 5th- order KdV equation, has been 
derived in (20] (with the parameter M defined in a different way) as a model equation 
for the gravity-capillary shallow water waves of small amplitude when the Bond number 
is close to but just less than 1/3. It has been extensively studied since then, see e.g. 
[21], and, although it is not integrable via the inverse scattering transform, it is known 
to have a rich structure of solitary wave solutions - in particular, existence of nonlocal 
solitary waves with propagating oscillatory tails and of asymmetric solitary waves has 
been established. Our analysis shows that the 5th-order KdV equation ( 1331) arises as 
the leading order equation in the classical water wave problem without surface tension 
when a = 0{(5 2 ). 

We will also present without derivation the leading order equation for the case 
= 0(a 3 ) obtained by retaining the terms which are at most cubic in a. It reads 

m + Vx + ^ar]r] x - ^a 2 r] 2 r] x + ^a 3 r] 3 r] x + ^tj 3x = 0. (34) 
This equation can be transformed into 

Vt = V 3 Vx + Vsx, (35) 

which belongs to the type K(m, n) introduced by Rosenau and Hyman [22J with m = 4 
and n = 1. Equation ( 135]) is nonintegrable but admits soliton-like traveling wave 
solutions in some range of wave velocities. 



5. Discussion 



We have presented a procedure for systematic derivation of the leading order and 
higher order evolution equations for the surface elevation of unidirectional shallow 
water waves. This procedure is based on a consistent ordering of terms in the original 
asymptotic expansions for a prescribed relationship between orders of magnitude of two 
small parameters a and 0. Our results provide a set of consistent model equations 
for unidirectional water waves which replace the KdV equation and the higher order 
KdV equations in the cases when the parameters a and (3 are not of the same order 
of magnitude. Some of the equations emerging in our analysis as the leading order 
equations in the asymptotic expansion for the unidirectional water waves have been 
proposed before as model equations in other physical contexts (e.g., the Gardner 
equation, the modified KdV equation, and the so-called 5th-order KdV equation). In 
the higher orders of approximation, a variety of evolution equations which can serve 
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as higher order models for unidirectional water waves on equal footing with the higher 
order KdV equations are found. Our analysis also reveals that certain model equations 
used in the literature are questionable since they have been obtained as a result of an 
improper ordering which is invalid for any relationship among orders of a and /3. 
The present analysis is based on assuming a prescribed relationship between orders of 
magnitude of two small parameters a and /3. However, the results can be interpreted 
in another, alternative way, along the lines of the analysis presented in [23J. The main 
concern of the analysis of |23j is to demonstrate that the condition /3 = 0(a) is not 
necessary for having a balance between nonlinearity and dispersion characteristic of the 
KdV equation and that the KdV balance is possible for any (3 provided that a — > 0. To 
this end the variables are transformed in such a way that the parameter /3 is scaled out, in 
favour of a, which leads to a prescription, in asymptotic terms, of the region of time and 
space where the balance occurs and so the KdV equation is valid. This conceptual shift 
from a relationship between orders of magnitude of the two small parameters to distances 
and times needed for achieving the balance between nonlinearity and dispersion provides 
a new view which is more relevant to applications in nature. 

However, the analysis of [23] is restrictive in the sense that the transformation 
of variables introduced in [23] may result only in the problem which leads to the 
KdV equation to leading order as a — > 0. In what follows, we show that it is not 
because of some intrinsic properties of the water wave equations but simply due to 
a specific character of the transformation used in [23J. We extend the analysis of 
[23] by introducing a generalized transformation dependent on a parameter n (the 
transformation of |23j becomes a particular case). This generalized transformation, 
like the transformation introduced in [23], results in the system of equations which 
contains only one small parameter a. Specifying the transformation parameter n to 
different values allows to obtain a variety of different problems and a variety of the 
corresponding leading order equations (like the Gardner equation, the 5th-order KdV 
equation and so on) including the KdV equation. As a matter of fact, each problem 
obtained from the original one by applying the transformation for a particular value 
of n can be equivalently obtained by assuming the relationship /3 = 0(a n ) between 
orders of magnitude of the small parameters. The former approach allows to consider 
the nonlinear-dispersion balance, epitomized by the soliton equations, as existing for 
any (3, provided that a — > 0, but imposes conditions on the regions of space and time in 
which the soliton dynamics (the KdV dynamics, the Gardner dynamics, the 5th-order 
KdV dynamics and so on) are expected to occur. 

In [23], the transformations eliminating j3 are applied to the original system of 
equations in terms of velocities (u,w), pressure p and elevation r] and then the system 
of equations in terms of and r] is obtained from the transformed equations. Therefore, 
in our analysis, we will also deal with the original equations (although the same could 
be done for equations (E]) (J5]) in terms of and rf). The system of equations of a 
two-dimensional irrotational wave, with effects of surface tension negligible, after non- 



Ordering of parameters in the water wave problem 



12 



dimensionalising takes the form 

u t + a {uu x + wUy) = —p x , (3 (w t + a (uw x + ww y )) = —p y , (36) 

u x + w y = 0, u y - (3w x = 0, (37) 

w = on y = 0, (38) 

p — r), w = Tj t + aurj x on y — 1 + an. (39) 

(In the notation of [23], y — >■ z, (3 — > 5 2 , a — > e.) The scales for (x,y,t) are as in (J6]) 
and the scales for u, w and p are respectively (a/ H)^gH, (a/ ' L)^/gH and pga. 
The following transformations are applied to equations f )36|) - f )39|) in 



a/5 \/H -v/a 

x ^ — x, y^y, t-t-=t, p^p, v ->Vi w -> ^= w (40) 

Va Va VP 



As the result the system (l36j)-(139j) reduces to the system of equations 

u t + a [uu x + WUy) = —p x , a (w t + a {uw x + ww y )) = —p y , (41) 

u x + w y = 0, u y - aw x = 0, (42) 

w = on y = 0, (43) 

p — r], w — r]t + aM77 x on y = 1 + ar/ (44) 

which are the same as (I36 i) -( l39j) . but with (3 replaced by a, for arbitrary (3. From the 
analysis made in section 3 of the present paper (and from an equivalent analysis of [23J) 
it is evident that equations (l4Tl) -(l44l) constitute the representation that leads to the 
KdV equation f|T7|) to leading order as a — > 0. 

As it was explained above, the transformations ( |40|) can be generalized. The 
generalized transformations are 

x ~^ ~^j2 x ' y->v> 1 p p> v^v, w -> ^ (45) 

where n is arbitrary. Applying the transformations (|45|) to equations (!36l) - (l39l) results 
in the system 

u t + a {uu x + WUy) = -p x , a n (w t + a (uw x + ww y )) = -p y , (46) 
u x + w y = 0, u y — a n w x = 0, (47) 
w = on y = 0, (48) 
p = r], w = T] t + auTj x on y = 1 + ar] (49) 

which is (|36l)-(l39l). with (3 replaced by a". 

It is clear that the same problem (|46l) - (l49l) results also from the assumption 
(3 = 0(a n ) which allows to replace (3 by a n without loss of generality. As a matter of fact, 
all the results of the present paper are related to the problem fr46l) - fj49l . independently 
of the approach through which it is obtained. The analysis made in Sections 3 and 4 
indicates that equations (146l)-(l49l) specified to different values of n give rise to different 
equations to leading order as a — > 0, The KdV equation arises as a particular case for 
n — 1. Other particular cases might be the Gardner equation (n = 2), the 5th-order 
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KdV equation (n = 1/2), the K(4, l)-type equation in the sense of [22] (n = 3). In 
the case, when the system fj46l) - fH9l) is treated as obtained via the transformations fj45|) . 
the results are valid under some conditions on the regions of space and time where thus 
the corresponding soliton dynamics are expected to occur. It should be emphasized, 
however, that, although the system (H6l)-(j49l) can be equivalently obtained either by 
applying the transformations (j4"5l) or by assuming (3 = 0(a n ), these two approaches 
represent alternative views which cannot be combined. 
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Appendix. Results for the case of nonzero surface tension 

Since the multiplier (1 — 3r) appears in the coefficients of the highest derivatives in the 
leading order equations, the case of r = 1/3 should be considered separately. We will 
assume that r ^ 1/3 and moreover, that \r— 1/3| is not small. Indeed, if \r— 1/3| <C 1, 
one has to introduce yet another small parameter e = r— 1/3 and consider the asymptotic 
expansion with respect to e as well, see e.g. [20] . 



We will consider the case of = 0(a 2 ) keeping the terms that are at most quartic in 
a. Following the procedure described above, we obtain 



5.1. (3 = 0(a 2 ) 




arf a 2 rf f3 5 o 4 

W = T) 1 1 (2 — OT)ri2x OL 77 

' 4 8 6 V ; ' 64 ' 



) 



(50) 




360 



-19 + 30r + 45r 2 )77 5x = 



(51) 
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where z — J rf x dx. 

If we keep in (151]) the terms of order not greater than 0(a 2 ) to retain the dispersion 

and nonlinearity at the leading order then it reads 

3 3 B 

Vt + Vx + ^VVx ~ v<y 2 V 2 ' r lx + -r(l - 3r)r] 3x = 0. (52) 
z o b 



This is nothing but the Gardner equation which for r = coincides with equation (J2j 
discussed in Section 4. This equation can be reduced to the modified Korteweg-de Vries 
equation 

r]t = Kr] 3x + 6i] 2 r] x , 

where K = sign(/3(l — 3r)). 
5.2. B = 0(a 3 ) 

Keeping the terms up to the order of a 5 we have 

otri 2 a 2 ri 3 5a 3 r? 4 8(2 — 3r) 

w = 7? - H ' - + — -n 2x 

' 4 8 64 6 ' 

7a 4 77 5 aB ( , 77^ 

+ + T ( (2 + t)w2 * + (3 + 7r) T 

21aV + oV / 3(?T _ ^ _ _ + _ ^ 



512 16 V 2 

3 3 2 2 3 o o /3(l-3r) 

Vt + Vx + ~ a VVx ~ -at v Vx + —01 v Vx H 7 ^ 

2 8 lb b 

~~ l^" 4 ^ + ll^ 23 + 15t ^V2x ~ 2(3r - 5)77773^) + 7^ a V^ 

+ ^((" 13r + 19 )vl + 2(23 - 5^7777^ + 2(r + 5)t/V) = 0, (54) 

where again z — J rf x dx. 

If we consider the leading order equation, i.e., restrict ourselves to the terms which 
are at most cubic in a, then fl53l) becomes 



Vt + Vx + T^aVVx ~ ^a 2 V 2 Vx + ^-a 3 r} 3 r} x + § (1 - 3r)r] 3x = 0. (55) 
2 8 lb b 

This equation can be further transformed into 

th = Mr] 3 r] x + 773*, (5b) 

where M = sign(/3(l — 3r)). Eq.(JSBJ belongs to the type K(m, n) introduced by Rosenau 
and Hyman [22] with m = 4 and n = 1 . It is nonintegrable but it is readily seen to admit 
soliton-like traveling wave solutions for M = 1 in a certain range of wave velocities. 
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8 2 . 1 r r 2 arf 

w — 77 H — (2 — 3r)r?2x + P ( )??4:e 

' 6 V nx H K 10 6 8 ; ' 4 

(2 
(3 + 7r)^ + (2 + r)^£ 
lb 4 

4 / 1261 _ r _ r 2 _ _ 5r*\ 
V 113400 3780 80 48 128 J 

1 2^VV4x 



+ a/3 2 ( (326 + 435r + 315r 2 )^L + (28 - 20r + 5r 
V v 7 720 



720 v ' 80 

1091 + 480r + 945r 2 )^) + ^f, (57) 



8 3 8 2 

Vt + Vx + tt(1 - 3r)r7 3;c + -arr^ + —^-(19 - 30r - 45t 2 )t7 5x 
6 2 360 

+ |f,( 2 3 + 1 ^ + 2(5 -3^, 

' -275 + 399r + 315r 2 + 945t 3 )t) 7x 



15120 

_ ^a 2 r) 2 r] x + a/3 2 ((317 + 270r + 441r 2 )^|^ + (1079 - 150r + 855r 2 ) 
-(-57 + 50t + 15t 2 )||) 

+ ^ ^ - = 0. (58) 

P V 1814400 6048 2880 96 128/ ' v ; 

Notice that in the second order in 8, which is the leading order in this case, the 

equation (158j) for r] can be transformed into the Korteweg-de Vries equation only if 

19 - 30r - 45r 2 = (r w 0.4) when the term with r] 5x vanishes. If 19 - 30r - 45r 2 ^ 0, 

then ( 1581) in second order in 8 reads 

+ + 3r)r} 3x + \ocrm x + ^-(19 - 30r - A5r 2 ) V5x = 0. 

6 2 360 

We can get rid of the term r\ x by passing from x to x' = x — t, so upon omitting the 

prime at x the equation under study becomes 

Vt + ~(1 - 3 r)?fe + ^arjrj x + J- (19 - 30r - 45r 2 )r/ 5x = 0. 
6 2 360 

Next, let x = Ax, t = Bt, n = Cfj, where B = 360A 5 /(/5 2 (— 19 + 30r + 45r 2 )), 
C = -/3 2 (-19 + 30r + 45r 2 )/(1080aA 4 ). Then upon omitting tildes at x,t,r/ we obtain 
the so-called 5th-order KdV equation (see Section 4) in the form 

Vt = VVx + Kr] 3x + r] 5x , (59) 

where K = (3r — l)60A 2 /(/3(-19 + 30r + 45r 2 )). Assuming r ^ 1/3, we can set 
A = {\{8{-19 + 30r + 45r 2 )/(60(3r - 1))|) 1/2 , and then K = sign((/3(-19 + 30r + 
45r 2 )/(3r-l)). 
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5.4. a = 0(fB 3 ) 

If a = 0(/3 3 ), we obtain 



w = 7] + ^(2- 3r)r] 2x + (B 2 ( ^- - T - - T — ) rj 4x 



arf B 3 
~T~ ~ 15120 



(-488 + 756r + 630r 2 + 945r 3 ) 7] 6x 



nA ( 1261 61r r 2 r 3 5r 4 \ a/9 „ „ . 2 . 

+ ^ TToTnn - ^ " ™ " 7^ " TTTo )V8x + t^((3 + 7t)t7 2 + 4(2 + r) m2x ) 



V 113400 3780 80 48 128/ ' 16 

-159264 + 221936r + 161040r 2 + 249480r 



39916800 

+ 519750r 4 + 1091475r 5 )r7 102: + (l8(28 - 20r + 5T 2 )r]r] 4x 

+ (1091 + 480r + M5r 2 )r] x rj 3x + 2(326 + 435r + 315^)^) , (60) 
rh + Vx + ^(l- 3T) mx + |^(19 - 30r - A5r 2 )r] 5x 
+ \ot mx - (-275 + 399r + 315r 2 + 945t 3 )t7 7;e 



11813 55 19 2 r 3 5r 4 \ 



JV9x 



+ B 4 \ r t*- 

' V 1814400 6048 2880 96 128 

aB 

+ ^ (2(3r - 5)77773, + (23 + 15r)r] x r] 2x ) 

+ (95265 - 129943r - 90750r 2 - 131670r 3 - 259875r 4 - 1091475r 5 ) m lx 

39916800 v ; ' 



+ a/3 2 ((317 + 270r + 441r 2 ) ^M^L + (1079 - 150r + 855r 
V 288 



2\ VxV^ol 



1440 

-(-57 + 50r + 15r 2 )||). (61) 

In this case the leading order is three, and upon introducing a new variable x' = x — t 
the leading order equation for 77 reads 

Vt + |(1 - 3t) V3x + |^(19 - 30r - A5t 2 ) V5x 
3a B 3 

+ - 7j^(- 275 + 399r + 315r ' + 945r3 )^ = 0. (62) 
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